Exact Jastrow-Slater wave function for the one-dimensional Luttinger model 
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We show that it is possible to describe the ground state of the Luttinger model in terms of a 
Jastrow-Slater wave function. Moreover, our findings reveal that one-particle excitations and their 
corresponding dynamics can be faithfully represented only when a Jastrow factor of a similar form is 
applied to a coherent superposition of many Slater determinants. We discuss the possible relevance 
of this approach for the theoretical description of photoemission spectra in higher dimensionality, 
where the present wave function can be straightforwardly generalized and can be used as a variational 
ansatz, that is exact for the ID Luttinger model. 
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I. INTRODUCTION 

Recently, much progress has been done for under- 
standing the crucial role of strong electron correlation 
in photoemission spectra, namely the properties of the 
one-particle excitations and the corresponding dynamical 
Green's function*^ This topic has attracted a renewed 
attention due to the impressive progress in the energy 
and momentum resolution of angle resolved photoemis- 
sion experiments, which confirm spin charge separation 
in quasi one-dimensional systems^, one of the most im- 
portant effects induced by strong electron correlation. 

From the theoretical point of view, the pioneer work 
for the dynamical Green's function dates to 1970^, which 
showed the absence of coherence in the strongly cor- 
related regime of the Hubbard model. Later the self- 
consistent Born-approximation 5 - was introduced, provid- 
ing a surprisingly accurate description of the single hole 
dynamics in lattice models, such as the t-J model rel- 
evant for High-Tc superconductivity. In this context, 
P.W. Anderson suggested that the Fermi liquid picture 
could be violated not only in one dimension but also in 
higher dimensionality (D), and especially in 2D. From 
this speculation, a tremendous amount of work has been 
devoted to the subject, starting from the reconsidera- 
tion of the dynamical properties of a single hole in a 
quantum antiferromagnel£ii£, and the detailed analysis 
of the dynamical properties of a single minority spin elec- 
tron in a bath of fully polarized electrons with opposite 
S pin£±°.aLi2 

Finally, an important progress was achieved by the Dy- 
namical Mean Field Theory (DMFT), that was able to 
describe the important feature of the Kondo resonance 
when approaching metal insulator transition of the half- 
filled Hubbard model in infinite dimensions ^ 3 : 14 In this 
case it was also shown that the single-particle excita- 
tions in the proximity of the Mott transition may be 
highly nontrivial even in the metallic side. In particu- 
lar at low energy, a Kondo resonance appears between 
the expected upper and lower Hubbard bands and deter- 
mines the coherent quasi particle weight of the Fermi liq- 
uid metal, which vanishes exactly at the metal-insulator 



transition. The predictions of DMFT have been con- 
firmed by many experiments. For instance recently, Mott 
transition in Vanadium Oxide was clearly explained 15 . 
It represents also a theory capable of characterizing the 
two energy scales found in photoemission experiments of 
HTc compounds. 16 However, the dynamical properties 
predicted by this theory are not well understood out- 
side the DMFT formalism. In particular, it should be 
very important to characterize the anomalous low-energy 
excitations determining the Kondo resonance from the 
direct solution of the Schrodinger equation, namely by 
direct inspection of the cigenfunctions of strongly cor- 
related models such as the Hubbard model or the t-J 
model. 

In this work, we consider a much less ambitious task 
and we use an approach that can provide useful insights 
on the exact ground state wave function and excitations. 
We consider the well known Jastrow-Slater wave func- 
tion that has been used successfully in several corre- 
lated systems and we focus our analysis in one dimen- 
sion where analytic calculations are possible and numer- 
ical works have confirmed the impressive accuracy of the 
Jastrow-Slater wave function on several strongly corre- 
lated model a 17 : 18 ' 19 . This wave function can be generally 
written as a product of a Slater determinant, character- 
izing free electrons, times the so-called "Jastrow factor" 
J that appropriately weights the electron configurations 
(e.g. suppressing the wave function amplitude when the 
electrons are too close), in order to describe electron cor- 
relation. Indeed, by a lengthy but straightforward deriva- 
tion, we found that the exact ground state of the Lut- 
tinger model can be written as a Jastrow-Slater wave 
function. Moreover, not only the ground state, but also 
single-particle excitations can be written in a suitable 
Jastrow-Slater form. It is interesting that in this case, 
many Slater determinants have to be considered, with 
an appropriate change of the Jastrow factor. The results 
we obtained arc in perfect agreement with the Luttinger 
liquid theory in ID and the extension of the wave func- 
tion excitations to higher dimensions seem to imply the 
same Kondo-like resonance scenario obtained be means 
of DMFT, although at present this is just a speculative 



conclusion and a further numerical work is necessary to 
verify it. 

This paper is organized as follows: In section [ill we 
review the Luttinger model and later show how the one- 
dimensional fcrmionic Luttinger model Hamiltonian can 
be reduced to a quadratic Bose Hamiltonian by means of 
the so-called bosonization technique. In section IIII1 we 
show how the standard Bogoliubov transformation can 
be used to diagonalize the quadratic Bose Hamiltonian 
in order to obtain an exact ground state of the Luttinger 
model. In section HVl we show that the ground state of 
the Luttinger model can be rewritten as a Jastrow wave 
function, whereas in section [V] the single-particle exci- 
tations of the Luttinger model are given with an explicit 
Jastrow multi-Slater wave function. Finally in section 
Vll we study the dynamics of a single particle added to 
the right branch of the Fermi sea. 



II. FORMALISM 

We consider the one-dimensional Luttinger model fol- 
lowing the work of Lieb and Mattis^ with the notations 
given in a later work 21 : 



-v ¥ (k+k ? ) 




FIG. 1: Linearization of the realistic dispersion of a one- 
dimensional Fermi system about the Fermi points k F and 

— k F - 
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u 



(5) 
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where v F is the Fermi velocity describing a linear band 
around the Fermi momentum k F and u is a generic 
interaction which depends only on the distance be- 
tween electrons. In this linearization scheme (sec fig- 
ure [T|), the allowed momenta in the right (+) and left 
(-) branches satisfy the usual quantization conditions 
k = -j- x any integer , valid for periodic boundary condi- 
tions, assumed here and henceforth. These two branches 
are then extended to — oo, +oo by means of two fermionic 
fields ip± with their appropriate Fourier transforms 



After this important observation, it is possible to repre- 
sent these density operators in terms of canonically con- 
jugate bosonic fields <I> and X&.: 



N+(x) = -4- [Il(x) + d x $(x)) 



/4-K 



N-(x) = --U\n(x)-d x *(x) 



(6) 



with 



^± (&) = -7= / dx 



^±(x) 



(2) 



These fields define a local charge operator N(x) — 
N + (x) + iV_ (x) by means of the following contributions 
coming from the right and left branches: 



[*(a?), II(a;')] = iS(x - x') . (7) 

The total density operator N(x) can be written as 

N(x) = N+(x) + N_(x) = -^d x <S>{x) . (8) 

Now it takes just a little more algebra to show that the 
fcrmionic Luttinger Hamiltonian ([T]) can be expressed in 
terms of these bosonic fields where it assumes the follow- 
ing canonical harmonic form: 



N^(x) = ipl(x)np-(x) - (i>l(x)i/>-(x)) 



(3) 



h = yJ dxUn 2 (x) + ^[d x ^ 



(*)]' 



(9) 



The basic relations that make the Luttinger model ex- 
actly solvable are given by the nontrivial commutation 
rules of these left and right branches density operators 20 : 



W±(q), N ± (-q)] = T 



Lq 

2tt ' 



(4) 



where v s is the re-normalized Fermi velocity and K is 
the Luttinger parameter that can be simply expressed in 
terms of v F and the k = component of the interaction 
M i 20 i 23 . In principle, the mapping of Eq. (QJ to the har- 
monic Hamiltonian (J5|) is exact only if the interaction is 
assumed to be delta-like, but it can be easily generalized 
to a momentum-dependent coupling constant u, as it was 



done for instance in Ref^S. This canonical form given in 
Eq. (JHJ) embodies the entire low-energy physics of the so- 
called Luttinger liquids and due to its simplicity, can be 
solved explicitly. To this purpose, we introduce Fourier 
transforms of the bosonic fields: 

< 

n fe = -^ J L dxii(x)e ik * ~ u(x) = j- £ fc n k e-^ 

(10) 
By making use of the relations in fTO]) . ([9]) becomes: 



* = ?E 



k (n fe n_ fe ) + - (fc 2 $ fc $_ fc ) 



(11) 



which expresses the Hamiltonian in terms of "normal co- 
ordinates", $fe and life. Notice that we have essentially 
reduced the problem to a single harmonic oscillator for 
each given momentum. The next step merely repeats the 
procedure carried out for a single harmonic oscillator. We 
define a set of conjugate creation and annihilation oper- 
ators: 



*» 



IT 
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a\ + a-k 


y/2\k\ 


■ /W 
l \l 2 


al k - a k 



(12) 



Now substituting (|12[) into (fTTj) . we obtain the stan- 
dard quadratic Bose Hamiltonian: 



»=E 



Ck 



t 1 

a k a k + 2 



E? 



a!o_ fc + a k a- k 



where the functions 7/. and efc are given by: 



e* 



7fc 



#[*+*] 



_ v,\k\ 



[K-*] • 



(13) 



(14) 



The Hamiltonian (j 13() contains only excitations in a given 
sector of particle number because in this section, we are 
interested in determining only the ground state (GS). 



III. EXACT SOLUTION OF THE MODEL BY 
BOGOLIUBOV TRANSFORMATION 

We will use the following unitary transformation which 
is the standard Bogoliubov transformation to diagonalizc 
Eq. (USD: 



U = e ia , 
where S has been obtained in Ref. (|201): 

S = i 2^ 9q a q a -q + n - C - • 
9>0 



(15) 



It is easy to show that 

UdqW = cosh(9 q )a q + sinh{9 q )a_ 



(16) 



It follows immediately that the transformed Hamiltonian 
H — > UHW assumes the diagonal form 



UHrf = J2 ( E kaW 



k\E k >0 



<_k_ 

2 



where 

E k = y/el->yl=v a \k\ 

This can only be possible provided 

7fe 



tanh(6*fc) = — 



£fe + V4 - it 



2" • 



(17) 



(18) 



(19) 



namely that 9 k does not depend on the momentum k and 
is simply given in terms of the Luttinger parameter K, 



„-20>- 



K 



After this unitary transformation, the ground state 
cannot contain any boson excitation because E k > 0, 
implying that in this representation, the ground state co- 
incides with the vacuum |0), namely the non- interacting 
Fermi sea \FS), the unique state of the Luttinger model 
corresponding to the vacuum of the canonical operators 
a k (a k \FS)=0, Vfc). 

It follows that in the original representation, the GS 
of the Bose Hamiltonian is simply given by U\FS) and 
can be written in the following general form up to an 
irrelevant normalization constant: 






t t 

1 ' n ' 



\GS)=ex P nf q ay_A\FS) 

U>0 J 



(20) 



Here, f q is a function whose analytic form will be de- 
termined in what follows. In order to determine f q , we 
notice that UakU' must annihilate the GS U\FS) Vfc. 
After substituting U\FS) with the ansatz (f2"0"]) and by 
means of Eq. (|16|) , we obtain a simple equation 



(21) 



al fe (cosh(6 l fe)/ fc + sinh(9 k )) = 



which can be solved to get 

/ fe = -tanh(0 fc ) = 



Ik 



£fe + \ft 



Ik 



(22) 



By replacing the quantities tk and 7& from (jT4|) 
into (|22p , we obtain the pairing function in terms of the 
Luttinger interaction parameter K: 



fk 



1 



K 



(23) 



When K = 1 (i.e. in the non- interacting regime), the 
pairing function f k correctly vanishes, and the ground 
state is just the Fermi sea \FS). 



IV. GROUND STATE OF THE LUTTINGER 

MODEL BY THE JASTROW-SLATER WAVE 

FUNCTION 

In this section, we show that the ground state of the 
Luttingcr model is a Jastrow-Slater wave function i.e. 
the ground state (|20[) can be rewritten as a Jastrow wave 
function \ipj): 



\ijj)=e-^* v « N -> N -«\FS} 
where \FS) is the free Fermi sea, 



N„ 



z^$ 



7T 9 ~ Z<7 V 27r|gl v "" 



(al + a- q ) 



(24) 



(25) 



and v q are the momentum-dependent parameters. 

Henceforth we assume that v q > and that v q = 

L 

for q — 0, because the total charge J dx N(x) is con- 

o 
served. In this way we will obtain the analytic form of 

the Jastrow parameters (v q ) as functions of the pairing 
amplitude f q , by requiring that 

e-{E,>o««^«^-«}|i?S) =R a e&'»of ka t a -*}\FS} , 

(26) 
where R a is an overall constant which depends only on 
f q . In the following derivation we do not even need to 
assume that the pairing function /& is constant. 

Let us first introduce the Hubbard Stratonovich trans- 
formation for the Jastrow factor: 



N q -z* q N- q ) 



e -{E g >o^W 9 JV_ g } = f TT ^l e ~\z q \ 2 +^^(z q 
■* q>0 ^ 

(27) 
Using Eq. (|2l))) . the RHS of the above equation simplifies 
to 



J q>0 7F 



-\z q \ 2 +A q +B q 



where 



iq ]/^[^ Zq + a -i Z ^ 



B„ 



iq 



27r\q 



(a-gZq + dgZ*) 



(28) 

(29) 
(30) 



Now using the Baker Haussdorf Campbell formula 

e A q +B q = e A qe B qe -l/2[A q ,B q ] ( y&M if [A q ,[A q ,B q ]] = 

[Bq,[A q ,Bq]] — as in this case), the commutator in 
the previous expression can be explicitly evaluated and 
it is a constant: 



1 



^IM l- |2 



-[A q ,B q ]=- ] - 



(31) 



On the other hand, e Bq \FS) = \FS) because all non- 
vanishing powers of B q annihilate the vacuum. Thus, we 



obtain the following equation after applying the operator 
to \FS): 



(32) 



q>0 



where 



mv q ^ „ 



(33) 



By performing the remaining simple Hubbard- 
Statonovich transformation integral we obtain: 



IVo 



where 



R a e 



R a 



' ^g>0 14 



q>0 



,t„t 



•\FS) , 



o„ 



(34) 



Notice that this constant can be infinite if v q does not 
decay sufficiently fast for large q. The divergence in the 
infinite product can generally be removed by introducing 
a large momentum cutoff, e.g. |g| < A CM t, and taking into 
account that R a is just an overall normalization constant 
that does not change any physical expectation value even 
when the cutoff is sent to infinity. 

Now by a direct comparison with Eq. (|26p . we obtain 
that the Jastrow wave function ipj is the ground state of 
the Luttinger model if and only if 



<h, 



l + a a 



1 



K 



(35) 



It follows immediately from the above equation that for 
K < 1 , the ground state momentum-dependence of the 
Jastrow parameters can be expressed as 



■K I 1 



Remark 



(36) 



From Eq. (|3"6")) . it follows immediately that when K = 
1 (corresponding to the free theory), v q — V q and 
hence \ipj) reduces to \FS) which is the ground state of 
the free theory. 



V. SINGLE-PARTICLE EXCITATIONS BY A 
JASTROW MULTI-SLATER WAVE FUNCTION 

In general, not only the ground state but all the eigen- 
states of the Luttinger Hamiltonian can be written in a 



Jastrow multi-Slater form with appropriate Jastrow fac- 
tors. To this purpose, let us consider that an eigenstate 
of the non-interacting Luttinger model \FSk) (defined 
later) can be transformed onto an exact eigenstate of the 
interacting Luttinger model (K ^ 1) by means of the 
unitary transformation U — e tS (defined in the previous 
chapter in Eq. (fT5|) ). namely 



\<h 



e iS \FS k ) 



(37) 



Here, \FSk) is a suitable excited state in the free theory 
with an extra particle added to the right branch slightly 
above the Fermi momentum (as we are interested in low- 
energy excitations only). This excited state can be ex- 
pressed as 



\FS k 



dx 



-ikx n iA 



1>Ux)\FS) 



(38) 



The right-moving operator ^+(2^) after integration over 
x, creates a state, namely a non-interacting Slater deter- 
minant with an extra particle added to the right branch 
and with total momentum k + kp. 

The eigenstate \tpk) defined in (j3"7) can be expressed 
exactly in the following generalized Jastrow-Slater form 
(see appendix 1X1 for a detailed derivation): 



does not simply mean that one-particle excitations of the 
Luttinger liquid can be expressed as linear superpositions 
of Slater determinants (this implication would be trivial 
because any state can be expressed in this way) , but pro- 
vides an important restriction to the form of the wave 
function because the number of Slater determinants used 
in this ansatz remains much smaller than the dimension 
of the Hilbert space. 



B. Spectral property of the Luttinger liquid: 
quasiparticle weight 



the 



With this formalism, one can recover most of the ex- 
act results obtained with the conventional bosonization 
technique since we have represented the ground state and 
the one-particle excitations using a different (but equiv- 
alent) functional form. For instance, we evaluated the 
quasiparticle weight for the Luttinger liquid determined 
by the ground state and the lowest one-particle excita- 
tion. Our calculation obtained with Jastrow-Slater wave 
functions implies that the quasiparticle weight vanishes 
in the thermodynamic limit according to the power law 
Z ~ L- e , where = (K + K' 1 - 2)/2. This result tics 
in perfectly with Luttinger liquid theor y 20 ' 24 . 



\ipj,h) = I dx exp 




4x>,ivv_, 



(39) 



x exp 



Ev^, 



ipl(x)e- ,kx \FS) 



A. Remark 

In appendix \X\ we establish the dependence of h q on 
the interaction parameter K . It can be shown that when 
K = 1 (corresponding to the free theory), v g = 0, h q = 
and the excited state ipj t h reduces to \FSk) which is a 
single determinant. Generally, when K 7^ 1, h q 7^ and 
in this case the excited state ipJ,h represents a multi- 
determinant Jastrow correlated state, very similar to 
the Edward's ansatz used for the single spin flip state 
of the ferromagnetic Hubbard model described in the 
introduction 9 . Our derivation shows that this ansatz is 
exact for the Luttinger model. Strictly speaking, the 
number of determinants required to describe the ansatz 
is infinite because in Eq. (|39p . the variable x is contin- 
uous and a Slater determinant is required in general for 
each value of x. In practice however, we can go back to a 
lattice discretized version of the Luttinger model where 
the position of an electron is discretized and therefore 
for generic lattice models, the number of determinants 
required for describing the present ansatz should scale 
as the number of lattice sites. In this respect, our result 



VI. ONE-PARTICLE DYNAMICS 

The above result can be extended to represent the 
quantum mechanical time evolution e lHt of any one- 
particle state in a suitable time-dependent Jastrow-Slater 
form. 

In particular, we apply to the ground state written 

L 

as e lS \FS), the operator J dx e~ lkx tp\( x ) that creates a 


fermion with momentum k + kp in the right branch. In 

this section we determine the exact time evolved state 



!*(*,*)) 



JHt 



dx e- lkx ^ + {x)e lS \FS) 



(40) 



in terms of a generalized Jastrow-Slater form. This form 
is similar to what we obtained for the one-particle exci- 
tation \ipk) hi the previous section. Since all the single- 
particle excitations have the same two-body Jastrow po- 
tential v q , it is reasonable to expect the following ansatz 



|*Jhc(M)) - / dx C(t) e- ifcx exp 

" 



lJ2 v « N « N -« 



x exp 



J2h q (t)e iqx N q 



1>Ux)\FS) 



(41) 



Indeed this is the exact time-evolved state within 
the Luttinger model Hamiltonian provided the functions 



h q (t) are chosen appropriately 



see appendix |B|) . 



VII. DISCUSSION 

The most important outcome of this work is that it is 
possible to obtain an essentially exact description of the 
low-energy properties of one-dimensional correlated mod- 
els by means of a Jastrow-Slater wave function. Indeed, 
we have shown that the ground state of the Luttinger 
model; the well known and accepted model for describ- 
ing low-energy physics in one dimension, can be written 
exactly as a long-range Jastrow factor applied to the un- 
corrected Fermi sea. 

In addition, not only the ground state of the Lut- 
tinger model but also single-particle excitations can be 
described by the Jastrow wave function. In this case how- 
ever, the Jastrow factor is applied to many Slater deter- 
minants. Thus, unlike the non-degenerate ground state 
which is just a Jastrow-Slater single determinant, the ex- 
cited state has an intrinsic multi-determinant character. 
This multi-determinant state reflects the effect of insert- 
ing a single particle to the ground state of the Luttinger 
model. The wave function for the system with an extra 
particle changes drastically in form within the Jastrow- 
Slater ansatz i.e., from a single determinant to a multi- 
determinant state. The dynamical properties of the Lut- 
tinger liquid were also formulated within the Jastrow- 
Slater wave function. More specifically, the dynamics of 
a single particle added to the ground state of the Lut- 
tinger model can be expressed as a Jastrow multi-Slater 
state, with time-dependent and complex Jastrow factors. 

To avoid confusion, we remark here that this is just 
a new alternative point of view of the well known exact 
solution of the Luttinger model in one dimension. How- 
ever, it is important to emphasize that the Jastrow-Slater 
wave function approach can be easily extended to higher 
dimensionality and indeed, a long-range Jastrow factor 
applied to a Slater determinant has been used widely 
for electronic simulations based on the so-called quantum 
Monte Carlo technique 25 -. However, the present form of a 
multi-determinant Jastrow factor, that we have shown to 
be a necessary ingredient to deal with the single-particle 
excitation spectrum of the Luttinger model, has never 
been used before to the best of our knowledge. This 
form is particularly important in one dimension, in or- 
der to destroy the quasiparticle weight and determine the 
non-Fermi liquid behavior. Therefore, we expect that the 
extension of this wave function to higher dimensionality 
may lead to a deeper understanding of the photoemission 
spectrum of strongly correlated materials like high tem- 
perature superconductors which display unconventional 
behavio r 26 ' 27 . In these systems, the photoemission spec- 
trum is still controversial and unexplained. For instance, 
the strong momentum dependence in angle resolved pho- 
toemission experiment er 8 ' i 30 with Fermi arcs or hole 



pockets cannot be easily reproduced with the conven- 
tional single determinant Jastrow-Slater ansatz for the 
single particle excitations^ 6 - On the contrary, the exten- 
sion of the wave function defined in Eq. (|39|) to higher 
dimensionality provides a variational ansatz containing 
more variational freedom for the excitations (e.g. h q in 
Eq. 155]) . that may lead to more accurate results and 
possibly better agreement with experiments. Obviously, 
a systematic variational Monte Carlo study outside the 
scope of this work is necessary to confirm this interesting 
possibility. 

Finally, we would like to comment on the possible ex- 
planation of the Kondo resonance in the spectral weight 
of a metal, predicted by DMFT in infinite dimensions. 
In our approach, the wave function of an added particle 
with momentum k can be viewed as a coherent super- 
position of Slater determinants with a real space defect 
located at each space position x (see Eq. (I5TJ)) ). The 
excitation described in Eq. (|39[) is very similar to old 
types of wave functions^ and previous approaches to con- 
sider the single hole dynamics 2 -. The common feature of 
this approach with the Kondo problem is that the single- 
particle (or -hole) excitation acts like a real space impu- 
rity in the frame where the extra particle is taken fixed. 
The impurity problem is a peculiar characteristic of the 
Kondo model fixed point and therefore, we expect that 
the Jastrow-Slater approach should be able to introduce 
another energy scale, the Kondo one, in the problem of 
the photoemission spectrum, very similarly to the DMFT 
scenario. Again also in this case, numerical work is nec- 
essary to verify this issue. 
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APPENDIX A: DETAILED DERIVATION OF 
THE SINGLE-PARTICLE EXCITED STATE 

We report below the detailed derivation of the single- 
particle excited wave function. Although it may appear 
cumbersome and elaborated, it is indeed very simple con- 
ceptually. All one-particle excitations, as well as the 
ansatz state (I39|) . can always be normal ordered accord- 
ing to the obvious rule that: 

• the fermionic operator ip + (x) is always the leftmost 
one. 

• after that all bosonic terms can be ordered in the 
normal way: the creation operators a„ to the left 
and the destructions a q to the right positions. 

• In this way all the destruction operators a q disap- 
pear because they have to be applied to the vacuum 
and the final expression drastically simplifies. 

Indeed after the above three steps, it is easy to convince 
ourselves that one can generally obtain, both from the 
exact excitation Eq. (j3"T)) and the ansatz Eq. ([3"9")) . similar 
expressions containing, after the fermionic operator, a 
gaussian form of operators a* applied to the vacuum. In 
this way, it is possible to match the two expressions. 



1. First step, simplification of the ansatz \4>k) 

In order to apply e to the fermion field i/j + (x) we 
follow Ref.(l20f) that have derived that: 



where: 



e iS i>\{x)e iS = ip\_(x)exp 



where for q > 0, 



(Al) 



,LM 



-(cosh 9q — 1) 



■ sinh 



(A2) 



We can use Eq. (|Al[) to simplify the exact excitation \tpk) 
given by Eq. (|57jl : 



IV'fc) = / dx e- tkx e lS ^ + {x)e' lS e lS \FS) = (A3) 

o 

L 



J dx ip' + (x)exp 



ii 



Y,e iqx {h L q M a\ + h™a- q ) 



e lb \FS) 



We recall that the ground state of the Luttinger model, 
apart for a normalization constant, is given by: 



\GS)=e F \FS) 



(A4) 



t t 



q>0 



(A5) 



In this way we obtain: 

IV'fc) =fl«x 

L 



t („\ F C 



dxip + (x)e exp 



Finally by using 



£e^e-^4 + ^fa_ g y 



\FS) 



-F F _ f t i 

e~ F a p e F = a P 



(A6) 



and applying again e A+B — e A e B e 1 / 2 l A ' B ] with A = 
J2 q (hq M + fqh^y^al and B = J2 q h FM e~^ x a q we 
obtain: 



IV'fc) = C 



[ e~ lk2 



t <„\„F 



dx t\)\(x)e x 



x exp 



\FS) ,(A7) 



Y^{h L q M + J q h^)a\eJ 
C = c^(l/2Y,h L q M {h L q " + f q h^)\ . 

Second step: simplification of the ansatz \tpj 



Also the ansatz ipLj can be recast in a form similar 
to Eq. (|A7j) using even simpler algebra because the Jas- 
trow depends only on the total density operator N q , with 
commutation rules 



[N q ,N q ,] = 



N q ,^ + (x) 



vT c 



-*8*,/.t 



mx). 



(A8) 
(A9) 



In this way it is easy to derive the following useful rela- 
tions: 



-i/2J2 q v g N g N- q ^U x \ e i/2i:,«,iV,iV_, 



= ip\ (x) exp 



»*". + i^> 



e E, *.*.«'- ^ +{x) e -J2 q h q N q e^ 



— ip + (x) exp 



Ehq_ 

L q v 



With the above relations we can bring the operator 
tjj(x) — e lkpX, ipl r (x) (as all the left branch states are occu- 
pied for momenta p ~ kp as we assume) in the leftmost 
side of Eq. (|39|) and obtain: 



\i>jh) -- C Jh l dx ^ + {x)e~ ikx e F x 



\ipj,h) = exp 



E 



VI. 



x dx ip\{x)e lfea; exp 
o 



x exp 



Jewv- 



E4 2) ™, 



L q 



where 



and 



Z (D - h _ J!«_ 
9 g 2VI 



z (2) _ I, _ _^_ 
9 ? VI - 



(A10) 



(All) 



(A12) 



Our ansatz can be further simplified by implementing 
the condition which enabled us to obtain the Jastrow 
parameter in terms of the pairing function, namely : 



exp 



lJ2 v « N « N -« 



\FS)=R a e F \FS), (A13) 



where F has been previously defined in Eq. (|A5[) . Now 
we can replace N(q) in terms of canonical operators as 
in Eq. (|25[) and we can perform similar steps as before, 
namely: 



\ipjh) = i? Q exp 



Ji) 

L <? 



dx tp\_(x)e lkx x 



J 



x exp 



^ iqz q ] e lqx f 

V v 27r M 



e F |.FS) 



This can be also written as 



IV> 



Jh) 



i? Q exp 



x e^exp 



- Zq 

J 

zgzq ' e lqx 



dx ip\{x)e lkx x 



B5r> 

(2) 

Eiq^q c ■ _p, 4. . p 

— =^e (a\ + a_ q )e 

g V 27r kl 



\FS) 



x exp 



E 



a q 



V^M 



|F5)(A14) 



where 



C. 



Jh 



R a exp 



,(!) 



V^/^9 



l£l 
47r' 



(^ ili z (2) z ( 2 )Q I f ) 

^ \/L A ^ q q 

. q 



We are now in the position to match the two states 
\ij)jh) and l^fe), using equations (|A14|1 and (|A7|) respec- 
tively. Indeed, apart from an irrelevant constant, the 
ansatz \ipjh) is an exact excited state of the Luttingcr 
model, if the following condition is satisfied: 



iqz q 2 \l + f q ) 



= h L q M + f g h L _f. 



(A15) 



.(2) 



Since z q is linear in h q , the above equation is a sim- 
ple linear equation, that can be solved for the unknown 
quantity h q . 



APPENDIX B: DETAILED DERIVATION OF 

THE REAL TIME EVOLUTION FOR THE 

JASTROW SLATER WAVE FUNCTION 

1. Simplification of the state ^ jhc 

This can be obtained by applying the same derivation 
described in the subsection (|A2[) , with slightly different 
notations. Therefore the final expression is: 



\^j hC ) = I dxe- ikx i>l(x)C Jh (t)C(t) x 



x e F exp 



E 

L q 



v(2) 



igz q z >(l + f q )e* qx f 



V^M 



\FS) (Bl) 



C Jh {t) = R a exp 



y (i) 



^ ( VL 4tt q z - q[i + Jq) 



Here. 



and 



zfHt) = h q (t) - ^= 



(B2) 



Finally, using the relations in eq. (|A6|) and the fact that 
a q \FS) = Vq, we obtain: 



4 1} (*) = ^w ''" 



2y/L 



(B3) 



Simplification of the propagated state |\?t) 



where the constant M follows from the normal order of 



This is a quite cumbersome and complicated deriva- 
tion. We sketch how to obtain the final expression of 
this section. One writes the propagator: 

e iHt = e iS e iH t e -iS ( B4 ) 

Then one makes the effort to bring to the left the operator 
ift (x) using Eq. (jAip also for obtaining the expression of 
e~ lS ^ e lS which is the same result of Eq. (|Aip ( derived 
for e zS ip' + e~ zS ) with 8 q —> —0 qi namely 

h LM, 

Moreover we have to use that 



exp 



E 



e^ (h L q M a\ + h™a- q ) 



, the constant C\ follows from the normal ordering of 



\, \3q) -* h q ' = h q (~^«)- 



—F 

e exp 



E- 



,iq{x+v s t)f h LMJ , h LM 



(h q a q + h_ q a_ 9 ) 






1 



t ip{x— x')-\-iv s pt 



<,t r„i\ 



z lHot ^ + (x)e- lHot = / dx' [ ^J -e ip(x ~ x ' }+lVsPt I i>\{x') 

\ v ) 



i)\(x + v s t) 



e iHat a\e- iH ^ = e iE '*a\ 



and also the fact that 

e* s e* Hot ale- tHot e- lS = e im e iS a\e~ iS e~ im 



(B5) 
(B6) 

(B7) 



and C 2 (t) from the normal ordering of: 



— F 
e exp 



iq x h LM ( e «M cos/l (0 g ) a t + 



= e 1 ^ cosh{9 q )a\ + e^®* 1 sinh{9 q )a- q 

where E q = v s \q\ is not changed by the interaction. Then 
we arrive to the final result by using that e iH ° t \FS) = 
\FS), repeated applications of the Baker Haussdorf 
Campbell formula, and little extra effort such as 

e aa qe l3a\ _ ^a\ £ aa q g a/3 

(a relation that determines the constant C\ y2 {t) below): 

l*(M)) = 



MC X 



x exp 



C 2 (t)C l2 {t) I dx e- tkx ^ + {x + v s t) 


£ J"* ({h L q M + f q h^)e^"* + B q {t)) a\ 



L q 



x e \FS) where : 
B q (t) = h h q M [e iE ^cosh(e q ) + e- iE ^f q 8inh{e q )] 



Ci,2(t) = exp 

M = exp 

C\ — exp 

C 2 (t) = exp 



+e- ,;B '*sm/i(6l ? )a^ ? )] e i 



that can be made explicit by using simple manipulations 
already introduced in the previous section (see Eq |A6[) . 

It is clear therefore that, exactly as in the previous 
section, the above state can be written as a general- 
ized Jastrow Slater of the form \& jhCi with appropri- 
ate choice of the complex time dependent function h q (t) 
and time dependent constant C{t). Indeed after the 
simple replacement in the dummy integration in dx of 
x + v s t — > X (notice that we are using PBC and therefore 
iIj + {x + L) = ipAx)) so that the integration domain can 
be shifted by arbitrary constants), we obtain the follow- 
ing simple conditions to match: 



:lm\ 



^h^ M e~ lqVst B q {t) 
q 

1/2 5>j 

q 

l/2j2h L q M (h™ + h L _ff q ) 

q 

1/2 J2 eT lEqt h L J- q ' sinh{0 q )B q {i) 



C(t)C Jh (t) - e !tM MdC 2 (f)C 12 (t) 
iqz q 2) (t)(l + f q ) 



(B8) 



V^q] 



e - lqVs t ^ h LM + fqh LMyqv s t + g^ 

(B9) 



whereas the density density term in the Jastrow factor is 
always characterized by the same v q given by Eq. (|36p . 
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